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1.  Introduction 


The  electronic  properties  of  ferroelectrics  make  them  attractive  candidates  for  incorporation  into 
various  kinds  of  electronic  devices.  Because  the  performance  of  such  devices  in  a  given  circuit 
environment  depends  in  a  highly  specific  way  on  the  response  of  the  material  to  external 
perturbations  such  as  electric  fields,  stresses,  and  temperature,  it  is  important  to  provide  the 
device  designer  with  an  appropriate  set  of  parameters  that  characterize  this  response,  which  may 
or  may  not  be  easy  to  obtain.  In  contrast,  when  a  chemist  wishes  to  synthesize  a  new  material, 
he  needs  a  different  set  of  responses,  based  on  thennodynamics  and  chemistry.  Thus,  it  is  vital 
to  connect  these  differing  responses  mathematically  so  that  the  chemist  can  translate  data  from 
the  device  designer  into  a  form  that  allows  him  to  vary  the  parameters  in  the  proper  way. 

In  this  report  we  attempt  such  a  connection  for  the  specific  example  of  a  promising  new 
ferroelectric  material,  BaTi  9(Sc,Ta)  05O3.  This  material,  which  was  fabricated  by  a  novel 

process  developed  by  S.  Tidrow  and  the  first  author  of  this  report,  exhibits  a  strong  and  nonlinear 
ferroelectric  response,  namely  a  dependence  of  its  dielectric  constant  on  an  external  bias  voltage. 
Such  behavior,  which  is  referred  to  as  tunability  and  is  common  to  many  ferroelectric  materials, 
is  an  example  of  a  property  that  can  be  exploited  by  device  designers.  However,  because  it  is 
intimately  related  to  a  strong  increase  in  its  dielectric  constant  with  temperature  near  the 
transition  from  ferroelectric  to  paraelectric,  it  is  very  temperature-sensitive,  an  undesirable 
feature  for  electronic  device  applications.  Moreover,  most  procedures  that  decrease  the 
temperature  sensitivity  also  decrease  the  tunability,  so  that  it  is  important  to  determine  what  steps 
can  be  taken  to  decouple  these  processes. 


2.  Thermodynamics  of  Ferroelectrics 


The  analysis  used  in  this  report  is  based  on  the  macroscopic  theory  of  ferroelectricity  developed 
by  L.  D.  Landau  (1937)  and  A.  F.  Devonshire  (1949),  now  referred  to  as  mean-field  theory, 
which  postulates  the  existence  of  a  free  energy  per  unit  volume  of  ferroelectric  that  is  analytic  at 
the  Curie  temperature  Tq,  i.e.,  possesses  a  power-series  expansion  around  that  temperature: 

F  =  -a(T-Tr)D2  +-bD4  +-cD6 
2  \  U  4  6 

Here  D  =  so Eex  +  P  is  the  electric  displacement,  Eex  is  an  externally  applied  electric  field,  and  P 
is  the  polarization,  i.e.,  the  dipole  moment  per  unit  volume,  whose  nonzero  value  in  the  absence 
of  an  external  electric  field  defines  the  ferroelectric  state  (7). 
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In  a  “normal”  material,  the  minimum  value  of  this  free  energy  is  obtained  when  D  =  0.  The 
quantities  a,  b,  c,  and  Tq  are  referred  to  as  Landau-Devonshire  parameters.  They  are  all 

characteristic  material  properties,  which  the  material  scientist  can  vary  by  modifying  the 
synthesis  procedure,  and  are  assumed  to  be  temperature-independent  (2). 

In  the  presence  of  an  external  field,  the  displacement  can  be  computed  by  minimizing  this  free 
energy  with  respect  to  D,  yielding  the  polynomial  expression 

Eex=a(T-Tc^D  +  bDi  +cD5  (1) 


As  Eex  — »  0  this  polynomial  has  only  the  trivial  solution  D  =  P  =  0  when  T  >  Tq,  while  for  T  < 
Tq  it  has  several  nonzero  solutions  in  addition  to  D  =  0.  It  is  found  that  these  nonzero  solutions 
describe  states  with  a  lower  free  energy  than  the  zero  solution  when  T  <  Tq.  For  T  >  Tc  an 
additional  differentiation  yields  the  small-signal  dielectric  constant: 


s 
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It  is  this  quantity  that  is  usually  measured  in  e.g.,  capacitance  experiments.  For  small  values  of 
Eex  we  derive  the  following  expression  for  the  dielectric  constant: 


s 


T-T 


C 


which  is  field-independent  and  exhibits  a  characteristic  Curie-law  behavior  around  Tq. 


In  experiments  that  involve  moderate  electric  fields,  the  measured  values  of  s (Eex)  for 
paraelectrics  (T  >  Tc)  normally  can  be  fitted  to  the  following  trend  line:  if  s (Eex)  =  8qk (Eex), 
then 

K{Eex’T  )x<°’T)~A(T)Eex2  +  B(T)EJ  +C(T)Ej 


Relating  this  to  the  Landau-Devonshire  parameters  is  not  straightforward,  due  to  the  implicit 
nature  of  equation  1 .  To  deal  with  this  problem,  we  expand  the  displacement  D  in  powers  of  the 
external  field:  if  D  =  +  k^E^  +  k qE^  •  •  •},  setting  a(T—  Tq)  =  a  gives 
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E  —  ex  D  +  bD  +  cD 

=  asg  +k5^  +  k7^  •••|  +  6sq^  (k{E  +  k^E3  +  k^E3  +  KrjE^  •• 

+  cs05  |K^  +  k^E3  +  k^E3  +  k7 ^  ’ ' '  j 

=  asQK^  +  laSgK^  +  *s03k3)e3  +|ag^K^  +  3bsQ3K2K^  +  cSq3k3  jii3 
+  |asQKy  +36sq3k^k^  +3teQ3K^K^  +5c8q3k^k^  jii3  ••• 

The  last  term  can  be  neglected  since  it  is  not  needed  to  determine  a,  b,  and  c,  so 

/  3  3 \  3  /  37  55\5 

E  =  asQK^+  oiSqK^  +b&Q  \E  +  asQK^+3bsQ  k^k^+cSq  \E  +• 


Equating  the  coefficients  of  like  powers  gives 

1  =  as^K^ 

3  3 

0  =  ocSqK^  +bsQ 

3  2  5  5 

0  =  asQK^ +3bsQ  k^k^+cSq 


s0Ki=- 
u  1  a 


£0k3  “ 


£0k5  “ 


a(coki 


oc  4 

^  a 


3h^SQK^  snK~, +cfa~K,  ^  tl2 


0^3^L  (s0Klj  36^  -ac 


a 


a 


This  gives 


D  = ' E0Kl  Rx -*(e0Kl f  £«3  +(3*%'<1  -d(<=oKl f  V*  + ' 


and  hence  the  dielectric  constant 


£(E)  =  — =  £0Ki  i-3h(£oKi)  £2+5(3b2s0K1-cj(s0K1)  E4+- 


so  that 
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K(o,r)  =  K1 

so^(r)  =  3MsoKi)4 

so5(r)  =  5(362soKi-c)(soKi)6 


Note  that  even  if  the  Landau-Devonshire  parameters  are  temperature  independent,  the 
coefficients  A(T)  and  B(T)  will  vary  strongly  with  temperature  around  the  Curie  temperature  due 
to  the  presence  of  high  powers  of  Ki. 


The  coefficients  s(0,7),  A(T),  B(T),  and  C( T)  can  be  obtained  by  fitting  the  experimental  curves. 
A  least-squares  fit  for  the  Curie  temperature  of  the  form 


a  =  a{r  -Tg 


1 

s0K(0,r)’ 


where  a  and  Tq  are  fitting  parameters,  is  described  in  appendix.  Once  this  is  done,  it  is  easy  to 
derive  the  other  Landau-Devonshire  parameters: 

<2> 


5  LOW]  6+\A(Tf%  5h(0,r)]  7  (3) 


3.  Fitting  Results  for  BaTi#9(Sc,Ta)  Q503 


One  effective  way  to  modify  the  electrical  properties  of  perovskite  compounds  like  BaTi03  is  to 
dope  them  with  transition  metals.  We  have  found  that  this  leads  to  materials  that  are  both 
temperature  insensitive  and  highly  tunable,  a  circumstance  that  our  group  has  investigated  for 
several  years.  We  use  the  following  procedure  to  measure  the  dielectric  constant  versus 
temperature  and  field.  Contacts  are  deposited  by  E-beam  evaporation  on  circular  samples  of 
-10.5  mm  diameter  and  0.5  mm  thickness  to  fonn  parallel-plate  capacitors.  The  contacts  are 
composed  of  layers  of  250  A  Ti,  1500  A  Au,  3000  A  Ag,  1500  A  Au,  in  that  order.  The  small- 
signal  capacitance  is  measured  by  an  impedance  bridge  in  the  temperature  range  -55  <  T  < 

120  °C  at  5  values  of  bias  voltage  across  the  capacitors  from  0  to  500  V.  The  bias  voltages 
needed  for  the  measurement  are  provided  by  a  Bertan  205B  high-voltage  D.C.  power  supply. 
The  small-signal  capacitance  is  given  by  the  fonnula 

C  ~  £0*1  i^bias )  ft 
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Vr. 

where  A  is  the  contact  area,  h  is  the  sample  thickness,  and  ^^■ms  =  °^as  ,  the  electric  field 

created  within  the  sample  by  the  bias  voltage  Vjjias,  was  computed  for  each  sample  capacitor. 

This  fonnula  was  inverted  to  giver  the  dielectric  response  at  each  bias  field. 

It  is  a  striking  property  of  this  material  that,  despite  its  chemical  similarity  to  BaTiC>3,  it  is 

unambiguously  paraelectric  over  the  entire  temperature  measurement  range.  In  contrast,  pure 
BaTiC>3  is  ferroelectric  over  this  range.  A  typical  experimental  curve  of  the  dielectric  constant 

versus  field  looks  as  seen  figure  1 . 


The  symmetry  with  field  reflects  the  cubic  nature  of  the  crystal  symmetry,  which  eliminates  all 
odd  powers  of  the  electric  field  E.  Comparing  the  shape  of  this  curve  with  the  expression 

s  (£)  =  8j  jl  -  3bs3E2  +  5  |3h2s1  -  c J  s^4  +  •  •  •  J 

shows  that  it  curves  upward  at  very  high  fields,  indicating  that  3 b^e\  >  c  in  our  materials. 

As  stated  above,  our  fitting  procedure  starts  with  the  expression 

K(Eex’ThK  (°’T)-A(T)Eex2  +B(T)EJ  +C(T)eJ 

where  k(0,2)  is  taken  directly  from  the  zero-field  measurement  and  A(T),  B(T),  and  C(T)  are 
derived  by  least-squares  fitting  the  difference 

HEex’T)  =  K(Eex’T)-A°E 
=  ~a(t)e  2  +  b(t)e  4  +  c(t)e  6 

v  >  ex  v  '  ex  v  >  ex 
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Once  this  is  done,  the  other  Landau-Devonshire  parameters  are  computed  using  equations  2  and 
3.  Figure  2  shows  a  typical  curve  fit  of  the  experimental  points  (large  dots)  to  this  function  at  a 
specific  temperature. 


Figure  2.  Theoretical  fit  for  BaTi  9(Sc,Ta)  05O3  dielectric  constant  showing 
experimental  data  points. 


7  V  —  m 

The  fitted  values  of  a  and  Tq  for  this  material  are  4.2  x  10  -  and  -58  °C.  For  comparison, 

C  —  °K 

the  Curie  constant  C  =  1.8  x  105  °K  for  single-crystal  BaTi03  given  in  Lines  and  Glass  (3),  plus 

7  V  —  in 

the  transition  temperature  Tc  =  120  °C,  imply  that  a  =3.7  x  10  -  for  pure  BaTiCfi  (on  the 

C  —  °K 

ferroelectric  side  of  the  Curie  temperature).  While  the  value  of  a  for  our  material  is  reasonably 
close  to  this  number,  our  fitted  transition  temperature  for  the  (cubic -tetragonal)  transition  is 
remarkably  low.  Because  pure  BaTi03  undergoes  three  phase  transitions,  the  highest  (cubic- 


tetragonal)  at  around  120  °C  and  the  lowest  (orthorhombic-rhombohedral)  at  -60  °C,  it  is 
tempting  to  assume  that  in  our  material  the  two  high-temperature  phase  transitions  have  been 
disrupted,  leaving  only  the  lowest-temperature  one  to  form  the  observed  peak  in  s.  This  is  in  fact 
consistent  with  the  behavior  of  PbTi03-  Flow  ever,  x-ray  data  indicate  that  the  crystal  structure 


of  our  material  is  cubic  at  room  temperature,  suggesting  that  the  highest  Curie  temperature  has 
indeed  been  moved  down  by  this  large  amount. 
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Figure  3.  Fitted  temperature  variation  of  parameter  a  for  BaTi  9(Sc,Ta)  05O3. 

The  other  LD  parameters  b  and  c  are  strikingly  different  from  those  of  pure  BaTiC^. 
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c  -  parameter 


Temperature 


Figure  5.  Fitted  temperature  variation  of  L.-D.  parameter  c  for  BaTi  9(Sc,Ta)  05O3. 


Enumerating  these  differences,  we  have  that 


1 .  According  to  LD  theory,  the  constant  b  should  be  negative  for  ferroelectrics  with  first-order 
phase  transitions  like  BaTiC^.  For  our  material  this  number  is  positive. 


2.  LD  theory  predicts  that  b  and  c  should  be  temperature  independent,  whereas  both  are 
strongly  temperature  dependent  in  our  material,  although  there  is  a  range  of  temperatures 
where  they  are  almost  constant. 


3.  Both  b  and  c  are  much  larger  than  they  are  in  BaTi03-  Their  minimum  values  in  the  “flaf 

1 5  V  —  m3 

regions  of  these  curves  (around  room  temperature)  are  -2x10  - - —  and 

C3 


21  V  —  nx 

-8x10  - - —  respectively.  Contrast  these  numbers  with  the  values 

C5 

5  Q 

n  V-m  17  V-m 

b  =  6X101-3 - - —  and  c  =  6X101 ' - 7 —  reported  in  ( 4 ). 

C3  C5 


One  possible  explanation  for  these  anomalies  is  based  on  the  commonly  held  idea  that  first  order 
transitions  in  ferroelectrics  are  caused  by  electrostriction.  When  a  system  develops  a 
spontaneous  polarization,  the  associated  internal  field  causes  a  spontaneous  strain,  which  adds  a 

term  to  the  free  energy,  where  Q  is  the  electrostrictive  constant  and  7  is  a  combination 

of  elastic  constants.  Because  this  term  is  negative  and  proportional  to  D^,  it  reduces  b  in  the  free 
energy  expression  and  can  even  make  it  negative,  which  converts  the  phase  transition  from 
second  order  to  first  order.  Then  the  argument  can  be  made  that  in  BaTiC>3  the  constant  b  is 

normally  large,  but  is  canceled  out  by  this  term,  thereby  converting  a  second-order  phase 
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transition  into  the  (weak)  first-order  transition  observed  in  experiment.  We  tentatively  advance 
the  hypothesis  that  the  elements  we  have  added  to  the  BaTi03  somehow  eliminate  the 

spontaneous  strain,  converting  the  crystal  structure  from  tetragonal  to  cubic  and  “exposing”  a 
large  but  previously  hidden  positive  value  of  b.  Similar  arguments  can  be  made  for  c.  This 
hypothesis  will  be  a  subject  for  future  research. 
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Appendix.  Least-Squares  Fit  for  the  Curie  Temperature 


In  calculating  the  fitting  parameters  it  is  convenient  to  use  the  inverse  relation 


=  7  1  r  =  a(T-Tr) 
8(0,  T)  '  C) 


Then  the  error  function  is 


A{a’Tc)  =  ^[ai~a(Ti~TC 


where  the  summation  is  over  the  experimental  data  points.  Minimizing  this  function  with  respect 
to  a  and  Tq  gives 

^{aJc)  =  -^[TrTc}2[ara[TrTc)]^ 

J-a(»,7c)  =  zH2[ai-.(r.-rc)]  =  o 

U1C  i 

Set.  T.-Tc  =aE  Ti~Tc 

=>  j  1  1 

T.- T(j 


J.aiTi-TcJ.ara  £T2 -21^  +  NTC2 
i  i  V  i  i 


Xa.=al[r.-rc  : 
i  i 


HT.-NTC 


Let  N  be  the  number  of  data  points,  and 


A  =  — la. 
*i  1 


B  =  —TT. 
N  i  * 


c  =— sr.2 
1 


M=—  Za-T. 
N  i  11 


Then 
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M-ATc=a\C-2BTc+Tcz 


a  =  ■ 


B -T. 


C 


(B-rc)(M-Arc 


=  mc-ibtc+tc1 


BM  -  (M  +  AB)  Tc  +  ATC2  =  AC  -  2  AB  TQ  +  A T^ 
BM-(M  +  AB)TC  =  AC-2ABTC 
BM-AC  =  (M-AB)rc 
BM-AC 


=>  Tr  = 

c  M-AB 


a  =  ■ 


A  (M-AB) 


[BM-AC]  B  (M-AB) -(BM-AC) 


B~t  M-AB  ] 
A(M-AB)  _  M-AB 

-ab2  +  ac  c-b2 
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